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Abstract 

The influence of quadrupole and hexadecapole residual interactions on rota- 
tional bands is investigated in a single-j shell model. An exact shell-model di- 
agonalization of the quadrupole-plus-hexadecapole Hamiltonian demonstrates 
that the hexadecapole-hexadecapole interaction can sometimes produce a 
staggering of energy levels in the yrast sequence; however, long and regu- 
lar A/=2 sequences are not obtained. The shell-model results are discussed 
in terms of the intrinsic deformations extracted by means of the self-consistent 
Hartree-Fock method. The angular momentum dependence of intrinsic quad- 
rupole and hexadecapole moments and is investigated. 
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I. INTRODUCTION 



Recently, there has been considerable interest in the behavior of hexadecapole deforma- 
tions at high spins, motivated by the observation of unexpected regular variations in the 
spectra of superdeformed (SD) bands. Namely, in several SD bands, 149 Gd |J, 153 Dy [Q, 
194 Hg [|, 192 T1 g§, 148 Gd §, and 131 > 132 Ce §, a AI=2 staggering of the dynamical moment 
of inertia, J^ 2 \ has been observed. The effect is very weak and in some cases rather uncer- 
tain J7|. It manifests itself in a systematic shift of every other state in a rotational band. 
That is, the sequence of states with spins J, 1+4, 1+8, etc., is shifted down with respect to 
the sequence 1+2, 1+6, J+10, etc. In other words, the transition energies of those bands 
can be parametrized as: 

£ 7 (/) = £ 7 (/) + e (/)(-l)i (1) 

where E 1 (I), the average (reference) 7-ray energy, and the perturbation e(I) are smoothly 
varying functions of I. According to the data, the amplitude of the staggering, |e(J)|, is very 
small - of the order of 100 eV. 

Since these oscillations separate a band into two families in which spins differ by four 
units of angular momentum, it seems natural to attribute their origin to a coupling between 
the rotational motion and hexadecapole vibrations This scenario requires the presence 
of a small component with a four- fold rotational symmetry, C4, in the mean field of the 
nucleus. 

The effect of such a term has been investigated in phenomenological models assuming the 
C4 symmetry axis coinciding either with the symmetry axis (z-axis) of the quadrupole tensor 
jj~T| or with the rotation axis (x-axis) ]1,12,13|. These theoretical studies have assumed 



the presence of a non-axial hexadecapole term in the Hamiltonian of the rotating nucleus. 
However, microscopic calculations based on the shell correction method [|14]-[17] or on the 
self-consistent Hartree-Fock (HF) method ||T8| , [l9l predicted an extremely weak collectivity 



associated with the non-axial A=4 fields. 

In Ref. pOj , in a simple quadrupole-plus-hexadecapole model which does not impose a 



priori any intrinsic deformations, it has been demonstrated that the stagg ering in may 
occur in certain cases. The main objective of the present study, based on the same model, 
is two-fold. First, we analyze the fluctuations (staggering) in the yrast line and discuss this 
phenomenon in terms of quadrupole and hexadecapole interactions. Second, we analyze the 
angular momentum dependence of intrinsic hexadecapole moments, Q^, especially those 
with \i=2 and 4. To the best of our knowledge, no systematic study of this effect exists 
in literature. Also, experimentally, very little is known about the non-axial hexadecapole 
deformations, Q42 and Q 44 , and about their /-dependence. Some limited evidence for these 
exotic deformations exists from the alpha scattering studies |2l|j2^] around 168 Er. 



The model to test the subtle phenomena discussed above should be capable of taking 
into account the interplay between rotation and shape dynamics of a many-body system. 
In our study, we investigate a model which can be solved exactly, namely, the model of a 
single- j shell filled with an even number of identical nucleons interacting via the multipole 
forces. This model, in the quadrupole-quadrupole variant El, was widely exploited in the 
past to study nuclear collective effects associated with quadrupole degrees of freedom. In 
spite of its simplicity and a rather limited configuration space, this model is able to describe 



2 



a variety of collective phenomena including collective rotation. For instance, even for as 
few as four particles in a single- j shell, Mulhall and Sips |24 found rotational structures 



present in the spectra. Collective spectra for various particle numbers and different values 
of j were analyzed by Friedman and Kelson [25]. Baranger and Kumar [^,0 studied 



quadrupole deformability and the influence of pairing, and Arima |28[] discussed excited 
rotational structures. Recently, within the large- j shell model, the quadrupole collectivity 
was studied by Burzyhski and Dobaczewski 



The paper is organized as follows. The single- j shell model is briefly described in Sec. |H|. 
Section contains the shell-model results. In particular, high-spin fluctuations in the 
yrast line are discussed in terms of the interplay between the quadrupole and hexadecapole 
components of the residual interaction. The corresponding intrinsic deformations are studied 
in Sec. |V| within the self-consistent Hartree-Fock (HF) approach. Finally, conclusions are 
contained in Sec. [V|- 

II. THE MODEL 

The model studied in this paper is an extension of a single- j shell model of Refs. fl23| , |24|| . 
Its Hilbert space is that of a (large) single- j shell occupied by identical N nucleons (iV-even). 
The model Hamiltonian is rotationally invariant and takes the form: 

H = H s . p . + V QQ , (2) 

where H s . p . is the single-particle Hamiltonian, while 

Vqq = -X2Q2 ■ Q2 ~ X4Q4 ■ Qa (3) 

represents the quadrupole-quadrupole and hexadecapole-hexadecapole residual interaction. 
For the single- j shell, H S . P . reduces to a constant assumed to be zero in the following. The 
multipole operators appearing in Eq. (|]) can be expressed as 

Q\v = X!O m i m 'IV) fl iv, (4) 

mm' 

where a m = Ta m T~ l (T is the time reversal operator) and the dot symbols in Eq. @ denote 
the scalar product 

Qx ■ Qx = J2(-TQ^Q^- ( 5 ) 

The Hamiltonian (fJFJ) is fully described by the two constants X2 and \i representing the 
strengths of the quadrupole and hexadecapole interactions, respectively. The dimension of 
the corresponding Hilbert space is finite. Consequently, the exact solutions of this quantum 



many-body problem can be obtained by means of a straightforward diagonalization |29 
This has been done in three steps, as described in the following. First, the M-scheme basis 
was constructed. It corresponds to iV-particle wave functions \NMn), with M being the 
projection of the total angular momentum /, and n distinguishing between states with the 
same value of M. Since the model Hamiltonian is rotationally invariant, its energies do not 
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depend on M. Hence M was assumed to be equal to zero. Second, the /-scheme basis \NIk) 
(M=0) was obtained by numerically diagonalizing the operator I 2 (/ is the total angular 
momentum of the system) in the M-scheme basis. Here, the additional quantum number k 
distinguishes between states with the same angular momentum I. Finally, the Hamiltonian 
has been diagonalized in the /-scheme basis. 

For the half-filled shell there appears an additional quantum number related to the 
particle-hole symmetry associated with the unitary symmetry operator 

W= II « + (6) 

m=~j 

[An arbitrary but fixed order of factors in Eq. (^) is assumed.] 

Since IA 2 ={— l)-^ 1 / 2 , depending on the value of j, IA is hermitian or antihermitian. By 
employing the identities: 

^ a m = - a mU, Ua m = -a+U, (7) 

it is easy to show that the particle-number operator N transforms under U as 

U + NU = 2j + l-N, (8) 

i.e., U transforms a state with iV particles into a state with 2j+l— N particles. Moreover, 
because of the relation 

U + Qxfi = -(-IfQx,-, (A>0), (9) 

IA commutes with Hamiltonian Consequently, the spectra of N- and (2j+l— iV)-particle 
systems are exactly degenerate, while for the half-filled shell, N—j+1/2, the eigenstates of 
Hamiltonian (H) acquire a new quantum number u=±i N , an eigenvalue of the hi operator. 

Properties of Hamiltonian @) depend only on the ratio x^ = Xi/X2- The absolute value 
of X2 determines the energy scale of the model. Consequently, in the following, we put 
%2=1 MeV. The "realistic" value of Xi can De estimated by means of the self-consistency 



relation of Refs. [|30"|-|32|j . In the harmonic oscillator approximation for the average potential 
and assuming the nuclear radius to be /2=1.2t4 1 / 3 fm, the self-consistent ratio Xi is equal to 
165/176, i.e., it is very close to unity. Therefore, in our study the value of %4 wa s varied 
from (no hexadecapole interaction) to 1 (self-consistent limit). 

The intrinsic deformations were obtained by means of the standard cranking HF proce- 
dure. The self-consistent one-body density in the rotating frame, p 1 ^, was found by solving 
the cranking HF equations 

h"(jr),p>] = o, (io) 

where If is the self-consistent Routhian, h UJ =h — ul x , h=Tr(VQQp UJ ) is the HF Hamiltonian 
(including the exchange term), I x is the component of the total angular momentum along the 
x axis, and to is the rotational frequency determined from the angular momentum equation: 
Tr(p"4)=/. 
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The maximum angular momentum carried out by a system of TV particles in a single- j 
shell is 



m + l-iV) 

'max 2 ' V xx / 

e.g., it is largest for a half-filled shell [/ max =(2j + l) 2 /8]. For instance, for j=19/2 and N=1Q, 

^max = 50. 

Having found the self-consistent density matrix, the intrinsic multipole moments, 

Q\fi = Tr(^Q v ), (A = 2, 4) (12) 



were calculated as functions of I. [The multipole operators appearing in Eq. (|12D are 
those of Eq. (§)•] Since it is assumed that the self-consistent density matrix has three 
symmetry planes (i.e., D 2 h is a self-consistent symmetry), the odd-/x components of Q2 and 
Q 4 vanish, and Qx^=Qx-^- This is consistent with the standard definition of the intrinsic 
system defined in terms of the principal axis of quadrupole tensor. 

It is convenient to relate the moments (A=2, 4) to simpler deformation parameters 
that guarantee the unique mapping of the Q-surface on the parameter values. The two 
quadrupole moments, Q 2 o and Q22, are usually written in terms of the "polar" coordinates 
/3 2 and 72 : 



Q20 = cos 72, Q22 = -^/5 2 sin72. (13) 

Since in our case p w has three symmetry planes, the hexadecapole tensor can be parametrized 
by means of the three "spherical" coordinates, /3 4 , 74, and £ 4 pjf : 

Q m = /3 4 I y ^ cos 5 4 + sin <5 4 cos 74 
1 

Q42 = —^PiSmSiSmji, (14) 
v2 

Q44 = & I \/^cos5 4 - y^sin5 4 cos7 4 

where 0<$4<7r and — vr<74<7r. Since in definition @ the radial dependence of Q has been 
ignored (the radial matrix element is the same for all magnetic substates, i.e., it is a scaling 
factor), the overall scale of the deformation variables /3 2 and /3 4 is given by the magnitude 
of the quadrupole operators (|T2"|), and is not related to the scale of the usual Bohr's shape 
parameters [[33] . 



(z) 

The hexadecapole shape that is axial with respect to the z-axis corresponds to 74=74 =0 
and 64=5^, where cos 5[ z) = J7/12 (4 z) ~0.7017 or 40°20')- Analogously, the hexadecapole 
shape which is axial with respect to the x-axis corresponds to 7 4 =74 r ' ) =27r/3 and 5^=S[ Z \ 
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III. SHELL MODEL RESULTS 



To analyze the staggering effect in collective bands, several criteria have been introduced 
in the literature. One possibility was suggested in Ref. [0 where the staggering was discussed 
in terms of 

A 3 £ 7 (J) = - [£ 7 (7 - 3) - 3E 7 (I - 1) + 3£ 7 (7 + 1) - E^I + 3)] , (15) 
while another staggering filter, 

(16) 



E 7 {I) - ^ (AE y (I - 2) + 4£ 7 (7 + 2) - E y (I - 4) - E y (I + 4)) 



was introduced in Ref. ||. The usefulness of quantities ( |i~5|) and ( |i~6|) for visualizing irregu- 
larities in rotational bands has recently been discussed in Ref. where it was concluded 
that, in some cases, they can dramatically overemphasize perturbations in rotational levels 
(which occur, for instance, due to accidental degeneracies). In this context, we would like to 
point out that one can construct a set of quantities A k E^, related to various derivatives of _E 7 
with respect to I. Indeed, for a function f(x) defined on a discrete equidistant grid, fi=f(xi) 
( x i+i — x i = h), its derivatives can be approximated by the finite difference formulas: 

(£) 4 (/■«-/,) (17a) 



x=(xi+x i+1 )/2 

dx 2 ) h 2 



- 2f t + f l+1 ) (17b) 



(g) « ^ (3/, - Vi+i - fi-t + M (17c) 

\ / x=(xi+x i+ i)/2 

~ ^ (6/, - 4 f m - 4/^ + fa + f t+2 ) . (17d) 

V / X=Xi 

It is easy to see that filters fllSD and (|I6| ) can be expressed as 

A^»-2(^) md A' B »(^). , 18 , 

That is, A 3 /^ and A 4 .E 7 are related, respectively, to the fourth and fifth derivative of the 
total energy, E, with respect to I. By going to higher derivatives in defining the staggering 
filter, one effectively spreads out the perturbation due to configuration mixing at I=Iq over 
many states around I=Iq. Consequently, in order to avoid this artificial effect, in the present 
study we employ the simplest filter, 

A 2 £ 7 (J) = 1 [E,(I - 2) - 2E 7 (I) + E,(I + 2)] , (19) 
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as a measure of fluctuations fl37fl . According to Eq. ( |17a| ), A 2 E^ is the second derivative of 
E 7 with respect to I, or 

A9rn d 1 1 dj® , x 

A ^ K2 57^i = - 2 (^— • (20) 

where is the second moment of inertia. It is worth noting that for a "C4" spectrum, 
[Eq. ([I])] with E 7 linear in / (perfect rotor) and e(/)=e=const., all filters A k E^ (k—2, 3, 4) 
give results proportional to e. 

We now proceed by applying the staggering filter fll9|) to the results of the exact shell 
model diagonalization of the single- j shell Hamiltonian (^j). Figure [1] displays the even-/ 
yrast lines for iV=8 and 10 particles moving in the j'=19/2 shell, and interacting with the 
quadrupole and hexadecapole interactions with X4=0, 0.5, and 1. For the half-filled j=19/2 
shell (N=10) the particle- hole symmetry U is preserved and the open full dots indicate 
states having u—1 and u—— 1, respectively. The corresponding values of A 2 E 7 are presented 
in Fig. |. 

In both cases, N=8 and 10, the yrast line at X4 = and 0.5 can be understood in terms of 
smooth collective bands connected by large stretched (Ai=2) transition matrix elements of 
the quadrupole operator, Eq. (|]). In the limit of the strong hexadecapole interaction, X4 = l> 
the yrast line becomes highly perturbed and rather irregular; the quadrupole collectivity is 
lost, i.e., the yrast states are no longer connected by large quadrupole matrix elements. 

For the half-filled shell, the yrast line obtained in the X4=0 and X4=0.5 variants consists of 
states with alternating values of u. Thanks to relation (|9|), multipole operators only connect 
states with different w's. Consequently, for small and intermediate values of %4, the yrast line 
consists of long quadrupole band sequences, i.e., states connected by the Q 2 operator, while 
the stretched (A/=4) hexadecapole transitions are strictly forbidden. The situation changes 
dramatically for large relative strengths xa- Here, due to strong hexadecapole correlations 
favoring the AJ=4 coupling, the alternating-M pattern breaks down and the neighboring 
states often have the same values of u, i.e., they are not connected by a multipole operator. 
Hence, the yrast line cannot be considered as a rotational band. For the N=8 case, u is not 
conserved, and the u — > —u selection rule for the transition operator does not hold. 

By comparing Fig. [I] with Fig. §, one can see that the staggering filter (|T9|) singles out 
all kinds of irregularities occurring along the yrast line. For example, the characteristic 
fluctuation at 7=16 and X4 = (f° r both N=8 and 10) results from a small kink in the yrast 
line caused by a band crossing. Due to large band interaction, this irregularity is barely 
visible in Fig. [I]. 

According to Fig. ||], values of A 2 i? 7 increase with xa- F° r X4 = l they are an order of 
magnitude larger than for x 4 =0.5. Relatively long sequences of regularly staggered points 
are visible for %4 = 1 at large values of the angular momentum. One has to note, however, 
that the quantity A 2 E 1 represents the staggering along the yrast line, which in the %4=1 
case does not correspond to a sequence of rotational bands. For ^4=0.5, the magnitude of 
staggering is smaller and the staggering sequences are much shorter. 

As discussed above, the yrast lines calculated for iV=10, and %4 = and X4=0.5 do consist 
of two alternating A/=4 sequences with the same values of u. However, this fact is not 
reflected in the corresponding staggering pattern. On the other hand, the yrast band for 
X4=l does not have a form of an alternating-w structure, but the staggering pattern is 
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clearly seen. It is, therefore, obvious that the sole existence of two different representations 
of states with opposite values of (— l) 7 / 2 is not a sufficient condition for the appearance of 
the staggering pattern. In any case, the hexadecapole interactions considered in our study 
do not differentiate between the two particle-hole symmetries which occur for the half-filled 
shell. 

From the shell-model results alone, one is not able to understand the origin of the longer 
or shorter staggering patterns obtained in the calculations. In order to analyze the underly- 
ing physics, we proceed in the next section by investigating the mean-field properties of the 
single- j shell model. 



IV. MEAN-FIELD APPROXIMATION 

The cranking HF equation (|10"D was solved for the Hamiltonian (0,0) and for the same 
values of j, N, and \i as used in the shell- model study. At the largest values of u, the HF 
solutions correspond to the fully aligned states, i.e., they are given by the Slater determi- 
nants of fermions occupying single-particle states with the largest available projections of 
the angular momentum on the rotation axis. This is illustrated in Fig. |3| which displays the 
single-particle routhians (the eigenvalues of Routhian, h UJ \h , )=e^\i')) and the corresponding 
single-particle angular- momentum alignments i v . We see that for the large rotational fre- 
quencies, the routhians can be represented by straight lines with slopes given by constant 
alignments. (Recall that in the HF theory the relation i v =—de^/duj holds only approxi- 
mately because of the cu-dependence of the average HF field.) 

The total alignments {I^—Y^u^v (the average values of I x in the many body HF states) 
are shown in Fig. [| for the same parameters as used in Figs. [3] and 0. For iV=8 and %4 = or 
0.5, the limit of a full alignment 7 max =48 is reached. A conspicuous feature of the cranking 
HF results presented in Fig. f| is the appearance of gaps in the calculated yrast line. This can 
be attributed to the old problem in the description of band crossing in terms of the cranking 



model [38]. Namely, in the region of band crossing the lowest self-consistent solution jumps 
as a function of u from one continuous family of states to another. 

In cases shown in Fig. 0, a band crossing occurs for X4=0 and 0.5 (and for both values 
of N). On the other hand, for %4=1 there appear many consecutive band crossings related 
to configuration changes, and the yrast line is composed of short pieces of continuous HF 
solutions. This effect reflects the fact discussed in Sec. [Ill] that for large values of \a the 
yrast states of the exact solution cannot be grouped in bands connected by strong quadrupole 
transition matrix elements. It also explains the strong irregularities seen in the exact yrast 
lines calculated with large hexadecapole interactions. When these irregularities are analyzed 
with a staggering filter, one may obtain sequences of staggered points; however, in the single- 
j shell model studied here this does not seem to be connected with any collective effect caused 
by the hexadecapole interaction. On the other hand, even an isolated band crossing, as seen 
in the HF results for X4=0 and 0.5, is reflected by a kink in the exact shell-model yrast line, 
i.e., short deviations of A 2 i? 7 from zero neatly exposed by the staggering filter, cf. Fig. |2|. 

The HF energies are compared in Fig. |l] with the results of the exact shell-model diag- 
onalization. One can see that for %4 = the mean field approximation reproduces the exact 



results well [29]. Also for X4=0.5 the approximation is fair, although the exact results visi- 



bly deviate at low angular momenta from a parabolic, rotor-like behavior given by the HF 
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solutions. In particular, for both values of band crossings appear in the HF curves very 
close to the kinks seen in the exact yrast spectra. 

On the other hand, for %4=1 the HF approximation fails to reproduce the exact states 
at low angular momenta. Interestingly, both shell-model and HF predict the existence of 
yrast traps at 7^16; the multiple band crossing phenomenon is clearly correlated with the 
staggered shell-model energies. 

At very high angular momenta, when the full alignment limit is reached, the HF results 
approach the exact solution. Indeed, for I—I max , the shell-model wave function is represented 
by a single, fully aligned Slater determinant, a HF state. 

Based on the HF solutions, we may now discuss the nature of crossing bands. To this 
end, we show in Figs. |] and ^ the average values of quadrupole and hexadecapole moments 
[Eq. (|I2"D1, respectively. One can see that the deformation patterns for %4 = and 0.5 are 
rather similar, and these patterns undergo qualitative changes when % 4 increases to 1. Since 
our aim is to discuss collective rotation, in the following we concentrate on the case of low 
hexadecapole strengths. 

Figure ^ shows the values of quadrupole deformation @ 2 and 72 [Eq. ([TJ)] for the HF yrast 
lines at X4=0.5. It is to be noted that in our paper we consequently employ the Bohr-Hill- 
Wheeler convention of 72 deformation, not the standard Lund convention. Consequently, 
the collective prolate rotation takes place at 72=0° (the z-axis is the symmetry axis) and 
7 2 =-120° (the y-axis is the symmetry axis), while the collective oblate rotation corresponds 
to 72=60° (the y-axis is the symmetry axis) or 72=±180° (the z-axis is the symmetry axis). 
The non-collective rotation around the x-axis corresponds to 72=-60° (oblate shape) and 
72=120° (prolate shape). The regions (-60° < 72 < 120°) and (72 < -60° and 72 > 120°) 
are physically equivalent. 

As seen in Fig. |7|, the results for N=8 and 10 are very similar. Namely, the quadrupole 
deformation @2 decreases with increasing angular momentum. This decrease is almost con- 
tinuous across the crossing point at about 7=16. On the other hand, at the crossing point 
there is a sudden change of the 72 values from about —180° (collective oblate rotation to 
about —150° (maximally triaxial shape with the intermediate axis being the x axis). 

In the fully aligned state, J=/ max , the equilibrium shape of a system with N<j+l/2 
corresponds to 72=72 r ' )= ~60°. This is so because the eigenstates of I x cannot have nonzero 
average values of (^2,^=21 i-e., the x axis must be the symmetry axis, and, moreover, for 
N<j+l/2 the shape must be oblate. In this limit, values of quadrupole moments along the z 

axis are proportional to the Q20X values calculated along the x axis. Namely: Q20 = — Q20X, 



Q22 = -y -Q20X, and P 2 = |Q20x|- This is seen in Fig. |7| for N=8 at the high-angular- 
momentum end of the band. On the other hand, for the half-filled shell the aligned state 
has /5 2 =0, and this point is approached along the 72=— 150° line. 

Of course, this description in terms of the shape characteristics pertains only to the 
average values of the quadrupole moments, Q20 and Q22, which have well-known values for 
ellipsoidal shapes. In the following, we refer to the bands below and above the crossing as the 
collective oblate and triaxial bands, respectively, and the band crossing can be interpreted 
as a sudden transition caused by a change in the shape of the system. 

We may now discuss the evolution of the hexadecapole shape of the system within the 
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oblate and triaxial bands. The oblate bands do not exhibit any tangible non-axial hexa- 
decapole deformations, i.e., they have Q^^Qaa^Q- On the other hand, in the triaxial bands 
all three components of the hexadecapole tensor are nonzero. In both bands the hexadeca- 
pole deformations are fairly similar for %4 = an d 0.5, i.e., they weakly depend on the fact 
whether the hexadecapole interaction has or has not been taken into account. This suggests 
a simple interpretation of the hexadecapole moments as being induced by the quadrupole de- 
formations. Indeed, ellipsoidal shapes have nonvanishing hexadecapole moments which can 
be expressed through quadrupole moments or deformations. This leads to rather stringent 
conditions for the hexadecapole moments 



The assumption of the strong coupling of A=4 deformations to the intrinsic frame defined 
through quadrupole deformations requires that the hexadecapole tensor depends on products 
of quadrupole tensors. This leaves a freedom of many possible parametrizations [^1] defined 



through three scalar functions h 2 , h 3 , and /i 4 , depending on fi 2 an d cos37 2 . In the case that 
only hs is different from zero, one obtains: 



^ 4 = |/l3l v(il + i^ cos2372 )' (21a) 

~ cos 372 , , , 

cos o 4 = ± . , (21b) 

J I + cos 2 372 

COS74 = ± COS72, SU174 = ± sin 72 . (21c) 

[There are two possible sign choices in Eqs. ( |21b| ) and (|21c|) : either 74 = 72 or 74 = 72 ± n.] 
Expressions fl2"ip define shapes for which the hexadecapole moments follow the quad- 
rupole moments and do not constitute any independent dynamical deformations. (For the 
general parametrization for the hexadecapole tensor, see Ref. |41|| .) According to Eq. ( J21a ) 



the value of \h 3 \ can always be chosen in such a way that the strong-coupling and self- 
consistent values of /3 4 are equal, i.e., /3 4 =/? 4 . On the other hand, the strong-coupling values 
of angles <5 4 and 74 depend solely on 72. 

It turns out that the HF results obtained for low and intermediate values of xa can be 
understood in terms of the strong-coupling limit given by Eq. (|2T|). Figure |8] shows the 
calculated HF hexadecapole deformations f3±, £4, and 74 together with the strong-coupling 
values of <5 4 and 74 given by Eqs. ( |21b| ) and fl21cp , respectively. On can see that the assump- 
tion of a strong coupling accounts very well for the complicated /-dependence of <5 4 and 74 
presented in Fig. |||. Moreover one should notice that for the case of N=8 the values of £4 
and 74 are approaching the limit of the fully aligned state along x-axis where only the value 
of Qaox remains nonzero. 

It is interesting to note that in the iV=10 variant and 7=50 (termination point), hexa- 
decapole deformation (3^ vanishes for x 4 =0.5. Indeed, in the limit of angular momentum 
alignment, the expectation value (jm x \Q\^\jm x ) is proportional to c^ (f )(jm x X0\jm x ) , i.e., 

it depends only on m 2 , for even values of A [(jm x X0\jm x ) ps P\ f 2 ^) and for large values of 
j\. Hence the HF expectation values vanish in the fully aligned state of the half-filled 
shell (J°P x (x)dx = for even values of A). 

In the case of N=8 particles, the average values of can be expressed by the axial 
hexadecapole moment calculated along the x axis. Namely, the following relations hold: 
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Qao = ^Qaox, Qa2 = —^ViOQaox, Qaa = -y —Qaox, and f3 4 = \Q40xl Consequently, for 

7=48 the parameters <5 4 and 74 are approaching the values of 8± and 'y^, respectively. 

We conclude this section by presenting in Fig. the hexadecapole deformation (3^ for 
N=8 and 10 as a function of the coupling constant xa at 7=40, i.e., for the triaxial bands. 
On can see that for both particle numbers, /3 4 steadily increases with xa in the whole range 
of xa- That is, the hexadecapole phase transition that would manifest itself in a rapid local 
increase of f3± with xa, is not present in the single- 7 shell model. 



V. CONCLUSIONS 

In the present paper we studied the single-? shell model describing identical particles in- 
teracting via the quadrupole and hexadecapole residual interaction. The model was solved 
exactly, and the resulting yrast structures obtained for different strengths of the hexadeca- 
pole interaction were analyzed. The model was also solved within the Hartree-Fock approx- 
imation, where the results could be interpreted in terms of the standard shape variables. 
Results of the calculations were presented for the 7=19/2 shell and particle numbers iV=8 
and 10. Calculations have also been performed for other values of j and N, but the main 
conclusions can be drawn from the restricted set of results presented in the paper. 

The purpose of our study was two-fold. First, we analyzed the conjecture that the 
hexadecapole A=4 degree of freedom might be responsible for the AI=2 staggering effect 
in rotational bands. Second, we investigated the evolution of the hexadecapole shapes with 
angular momentum. 

Our study shows that by including the hexadecapole interaction in the single-? model, one 
can sometimes obtain a staggering of yrast energies. However, a long and regular sequence 
of AI=2 staggered energies cannot be obtained. This negative result suggests that the 
experimental data probably cannot be explained in terms of the coupling between rotation 
and hexadecapole vibrations. 

For a relatively weak strength of the hexadecapole interaction, one obtains regular collec- 
tive quadrupole bands which do not exhibit any staggering. One also systematically obtains 
the effect of a crossing between two different shape configurations: an oblate-shape collective 
band at low angular momenta and a triaxial band at high angular momenta. The irregu- 
larities which are present in the yrast line in the crossing region should not be confused 
with the staggering phenomenon. As discussed in Ref. [3f|, such a misinterpretation may 



easily happen especially when a multi-point staggering filter is used to extract the staggering 
amplitude from the calculated energies. 

For a relatively strong hexadecapole interaction, the yrast line becomes rather irregular 
and the states are no longer connected by strong quadrupole transition matrix elements. 
In the mean-field picture, such a yrast structure can be interpreted as composed of several 
different bands crossing one another. In some cases discussed in this paper, such multiple 
crossings can give rise to staggering patterns. 

In the studied model, one obtains nonzero hexadecapole moments even without including 
the hexadecapole interaction. These moments are rather weakly affected by the hexadeca- 
pole interaction. The hexadecapole deformations calculated in the single-? model simply 



11 



follow the quadrupole deformations, i.e., they can be discussed as resulting from a strong 
coupling between the A=2 and A=4 modes. Within the collective oblate bands, the hexa- 
decapole moment is axial and oriented along the symmetry axis of the quadrupole deforma- 
tion. For triaxial bands, the hexadecapole deformation also becomes triaxial, i.e., all three 
components of the hexadecapole moment differ from zero. They can be nicely interpreted 
within the strong coupling assumption, and depend in a simple geometrical way on the 
quadrupole asymmetry angle 72. It is worth noting that our self-consistent results strongly 
prefer one particular parametrization of the hexadecapole tensor of Ref. |f£lf , namely the 
strong-coupling expression with h 2 =h 4: =0. 



ACKNOWLEDGMENTS 



This research was supported in part by the Polish Committee for Scientific Research 
under Contracts Nos. 2P03B 098 09 and 2 P03B 034 08, by a computational grant from the 
Interdisciplinary Centre for Mathematical and Computational Modeling (ICM) of Warsaw 
University, and by the U.S. Department of Energy (DOE) through Contract No. DE-FG02- 
96ER40963 with the University of Tennessee. Oak Ridge National Laboratory is managed 
for the U.S. DOE by Lockheed Martin Energy Research Corp. under Contract No. DE- 
AC05-96OR22464. 



12 



REFERENCES 



[1] S. Flibotte, et al, Phys. Rev. Lett. 71, 4299 (1993). 
[2] B. Cederwall, et al, Phys. Lett. B346, 244 (1995). 
[3] B. Cederwall, et al, Phys. Rev. Lett. 72, 3150 (1994). 
[4] S.M. Fischer et al, Phys. Rev. C53, 2126 (1996). 
[5] G. de Angelis et al, Phys. Rev. C53, 679 (1996). 
[6] A.T. Semple et al, Phys. Rev. Lett. 76, 3671 (1996). 

[7] R. Kriicken et al, Abstract, Int. Conf. on Nuclear Structure at the Limits, Argonne, 

Illinois, 1996, ANL/PHY-96/1, p. 21. 
[8] L.K. Peker, S. Pearlstein, J.O. Rasmussen, and J.H. Hamilton, Phys. Rev. Lett. 50, 

1749 (1983). 

[9] I. Hamamoto and B. Mottelson, Phys. Lett. B333, 294 (1994). 
[10] P. Magierski, K. Burzyhski, J. Dobaczewski, and W. Nazarewicz, Acta Phys. Pol. B26, 
291 (1995). 

[11] A.O. Macchiavelli et al, Phys. Rev. C51, Rl (1995). 

[12] I.M. Pavlichenkov and S. Flibotte, Phys. Rev. C51, R460 (1995). 

[13] P. Magierski, Acta Phys. Pol. B27, 127 (1996). 

[14] I. Ragnarsson, private communication, 1994. 

[15] S. Frauendorf, J. Meng, and J. Reif, Proc. Int. Conf. on Physics from Large 7-ray 

Detector Arrays, Berkeley, 1994, Volume 2, LBL-35687, p. 52. 
[16] F. Donau, S. Frauendorf, and J. Meng, Phys. Lett. B, in press. 

[17] W.D. Luo, A. Bouguettoucha, J. Dobaczewski, J. Dudek, and X. Li, Phys. Rev. C52, 
2989 (1995). 

[18] P. Magierski, P.-H. Heenen, and W. Nazarewicz, Phys. Rev. C51, R2880 (1995). 
[19] A. Bouguettoucha, J. Dobaczewski, J. Dudek, W. D. Luo, and H. Molique, to be pub- 
lished. 

[20] K. Burzyhski, P. Magierski, J. Dobaczewski, and W. Nazarewicz, Phys. Scr. T56, 228 
(1995). 

[21] I.M. Govil, H.W. Fulbright, D. Cline, E. Wesolowski, B. Kotlihski, A. Backlin, and K. 

Gridnev, Phys. Rev. C33, 793 (1986). 
[22] I.M. Govil, H.W. Fulbright, and D. Cline, Phys. Rev. C36, 1442 (1987). 
[23] S.T. Belyaev, Mat. Fys. Medd. Dan. Vid. Selsk. 31 (No. 11) (1959). 
[24] W.J. Mulhall and L. Sips, Nucl. Phys. 57, 565 (1964). 
[25] M.I. Friedman and I. Kelson, Nucl. Phys. A144, 209 (1970). 
[26] M. Baranger and K. Kumar, Nucl. Phys. 62, 113 (1965). 

[27] M. Baranger and K. Kumar, Nucl. Phys. A110, 490, 529 (1968); A122, 244 (1968). 
[28] A. Arima, Prog. Theor. Phys., Suppl. extra number, 489 (1968). 
[29] K. Burzyhski and J. Dobaczewski, Phys. Rev. C51, 1825 (1995). 
[30] D.J. Rowe, Phys. Rev. 162, 866 (1967). 

[31] A. Bohr and B.R. Mottelson, Nuclear Structure, vol. 2 (W.A. Benjamin, New York, 
1975). 

[32] H. Sakamoto and T. Kishimoto, Nucl. Phys. A501, 205 (1989). 

[33] A. Bohr, Mat. Fys. Medd. Dan. Vid. Selsk. 26, No. 14 (1952). 

[34] D.L. Hill and J.A. Wheeler, Phys. Rev. 89, 1102 (1953). 

[35] S.G. Rohozihski and A. Sobiczewski, Acta. Phys. Pol. B12, 1001 (1981). 



13 



[36] V. Reviol, H.-Q Jin, and L.L. Riedinger, Phys. Lett. B371, 19 (1996). 
[37] C. Baktash, private communication, 1995. 

[38] I. Hamamoto, Nucl. Phys. A271, 15 (1976); R.A. Sorensen, Rev. Mod. Phys. 45, 353 
(1973). 

[39] S.G. Rohozihski, Report INR 1520/VII/PH/B, Warszawa 1974, unpublished. 
[40] W. Nazarewicz and P. Rozmej, Nucl. Phys. A369, 396 (1981). 
[41] S.G. Rohozihski, to be published. 



14 



FIGURES 



FIG. 1. The exact yrast spectrum (circles) of the quadrupole-plus-hexadecapole model for the 
j=19/2 shell filled with N=8 (left) and 10 (right) particles, and X4=0 (top), 0.5 (middle), and 1 
(bottom). The yrast bands calculated in the HF method are indicated by solid lines. For iV=10, 
open and full circles correspond to the values u=l and u=—l, respectively, of the particle-hole 
symmetry quantum number. 

FIG. 2. Staggering parameter A 2 E^ [Eq. (|l9|)1 as a function of I along the yrast lines obtained 
in the shell-model results presented in Fig. |l|. 

FIG. 3. Top: Self-consistent single-particle routhians, e^, for j=19/2 and iV=10 (yrast line) 
as functions of the rotational frequency u). Bottom: the corresponding single-particle alignments. 
Solid lines: signature r = +i. Dashed lines: r = —i. 

FIG. 4. Angular momentum (the expectation value of I x in the cranking HF state) as a 
function of u for the j=19/2 shell filled with iV=8 (left) and 10 (right) particles, and X4=0 (top), 
0.5 (middle), and 1 (bottom). 

FIG. 5. Equilibrium quadrupole moments Q20 (solid line) and Q22 (dot-dashed line) for HF 
yrast solutions as functions of angular momentum for the j=19/2 shell with iV=8 (left) and 10 
(right) particles, and X4=0 (top), 0.5 (middle), and 1 (bottom). 

FIG. 6. Same as in Fig. |5| except for equilibrium hexadecapole moments Q40 (solid line) and 
Q42 (dot-dashed line), and Q44 (dotted line). 

FIG. 7. Equilibrium deformations /?2 and 72 [Eq. (|F])] as functions of angular momentum for 
the HF yrast states obtained in the j=19/2 shell filled with iV=8 (solid line) or iV=10 (dot-dashed 
line) particles, and for X4=0.5. The limit of axial quadrupole shapes are indicated by dotted lines 
(72^: axial symmetry with respect to x-axis; axial symmetry with respect to z-axis). 

FIG. 8. Equilibrium hexadecapole deformations P4, S4, and 74 [Eq. (|i~4|)1 as functions of angular 
momentum for j=19/2, X4=0-5, and iV=8 (left) and 10 (right). Values of 64 and 74 obtained in 
the strong-coupling limit [Eq. ([21])] are denoted by dashed lines. The limits of axial hexadecapole 
shapes are indicated by dotted lines (£4 , 7^'- axial symmetry with respect to x-axis; 5^ z \ 'y^: 
axial symmetry with respect to z-axis). 

FIG. 9. Equilibrium hexadecapole deformation (5^ as a function of X4 for the HF yrast states 
at a fixed angular momentum, 7=40. Calculations were performed for the j=19/2 shell filled with 
iV=8 (solid line) or iV=10 (dashed line) particles. 



15 




(A9IAI) A0H3N3 



